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We use a method of linearization to study the emergence of the future cosmological singularity 
characterized by finite value of the cosmological radius. We uncover such singularities that keep 
Hubble parameter finite while making all higher derivatives of the scale factor (starting out from the 
a) diverge as the cosmological singularity is approached. Since such singularities has been obtained 
before in the brane world model we name them the "brane-like" singularities. These singularities 
can occur during the expanding phase in usual Friedmann universe filled with both a self-acting, 
minimally coupled scalar field and a homogeneous tachyon field. We discover a new type of finite- 
time, future singularity which is different from type I-IV cosmological singularities in that it has the 
scale factor, pressure and density finite and nonzero. The generalization of w-singularity is obtained 
as well. 

PACS numbers: 98.80.Cq, 04.70.-s 
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I. INTRODUCTION 

Starting out from the discovery of the cosmic accelera- 
tion [l[ there have been constructed many models of the 
dark energy, including the very unusual ones: the phan- 
tom energy, the tachyon cosmologies, the brane worlds 
etc. Consideration of these models results in some un- 
expected conclusions about possibility of new cosmolog- 
ical doomsday scenarios: the Big Rip singularity (BRS) 
0, the Big Freeze singularity (BFS) 3, Q, the Sudden 
Future singularity (SFS) 5], the Big Boost singularity 
(BBtS) [6], and the Big Break singularity (BBS) 0, Q. 
In all these models the evolution ends with the curvature 
singularity, \ 'd(t) \ — > co, reachable in a finite proper time, 
say as t — > t s . BRS and BFS both take place in the phan- 
tom models but with the different equations of state. In 
particular, BRS takes place if w — p/p — const < — 1 
whereas BFS occurs for the dark energy in the form of 
a phantom generalized Chaplygin gas. Models with the 
SFS, BFS, BBtS and BBS singularities are characterized 
by a finite value of the cosmological radius but different 
values of Hubble expansion parameter H s = H(t s ) and 
different signs of (divergent) expression d s /a (cf. also 



00, 



a s < 00, 



a, < 00, 



-co, 



-co, 



H s = +00, — = +00, 
a. 



(BRS) 



(BFS) 



< H s < 00, — = -00, (SFS) 
a s 



a s < co, < H s < co, — = +00, (BBtS) 

(l s 



a s < co, H s = 0, — = -co. (BBS) 

a s 

Remark 1. One of classifications of singularities for the 
modified gravity was given in [Tol j ; for the classification 
and discussion concerned with avoiding the singularities 
in the alternative gravity dark energy models cf. [llj ]. 
Another classification of finite-time future singularities 
(Type I-IV singularities) is presented in [12J. According 
to this classification, the BRS is a singularity of type I, 
BFS is of type III, SFS and BBtS are type II and BBS 
- type II with p s = p(t s ) — (although this is a quite 
non-trivial special case of a type II singularities). Our 
classification doesn't contain singularities of IV type (for 
t — > t s , a —¥ a s , p — > 0, \p\ — > and higher derivatives of 
H diverge) but as we shall see in Sec. VI, the classifica- 
tion of Ref. [l2[ is not exactly complete too: the type IV 
is the special case of a more general type of singularities. 
Remark 2. Another type of "singularity" - so called w- 
singularity was obtained in [13}. This "singularity" has 
a finite scale factor, vanishing energy density and pres- 
sure, and the singular behavior manifesting itself only 
in a time-dependent barotropic index w(t). The w- 
singularities seem to be most similar to the type IV but 
are different nonetheless since they do not lead to any 
divergence of higher order derivatives of H [l3| . 

One surmises that ^-singularity is not a correct phys- 
ical singularity since all the physical values (i.e. den- 
sity, pressure and higher derivatives of the scale factor or 
Hubble roots) are finite. Moreover, the definition of w- 
singularity from the [l3| is an incomplete one. To show 
this let us consider the following form of the scale factor 



a(t) 



A(t s -ty 



(1) 
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([T]) is the special case of the general form of the scale 
factor from the [3 (with B = Q, A = a s , C/t n s = -A, 
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D = 1, n = m). One can show that for t — >■ t s : 

(a) type III singularity if < m < 1; 

(b) type II singularity if 1 < m < 2; 

(c) lo-singularity if m > 2. 

The case m = 1 results in a model with the constant 
barotropic index u> = —1/3. The case m — 2 is the most 
interesting one because 

4A 

p^O, p^ t^O, \w\ ->• oo. 

3a s 

and 

at t = t s . Thus we have some generalization of w- 
singularity such that the pressure is non-vanishing and 
finite at t = t s . 

The BRS and BFS have been obtained in the phan- 
tom cosmologies (BRS for the phantom perfect fluid with 
equation of state pj p — w = const < — 1 and BFS for the 
phantom Chaplygin models. Throughout the paper we'll 
stick to the metric units with 8irG/3 = c = 1). The BBtS 
is connected to the effect of the conformal anomaly that 
drives the expansion of the Universe to a maximal value 
of the Hubble constant, after which the solution becomes 
complex. The BBS takes place in tachyon models. 

Unlike BRS, BFS and BBtS altogether, the BBS and 
SFS are violating just the dominant energy condition 
(p > 0, — p < p < p). It is also possible to obtain 
some generalization of these singularities. In particular, 
generalization of the Sudden Future singularities (the so 
called Generalized Sudden Future singularities or GSFS) 
are singularities such that one has the derivative of pres- 
sure p(" l ~ 2 ) singularity which accompanies the blow-up 
of the m-th derivative of the scale factor These 
singularities are possible in theories with higher-order 
curvature quantum corrections [l2j and corresponds to 
classification in this paper. 

Despite the fact that there has recently been a great in- 
flow of articles, elaborating on the aforementioned singu- 
larities, an absolute majority of them has been of a math- 
ematical nature, while the physical reasons for arousal 
of such singularities still remain less then clear. A re- 
markable exception is the article [H, which has intro- 
duced for the first time a new type of cosmological sin- 
gularities located on the brane (for discussion about the 
soft singularities on brane with the quantum corrections 
cf. [HI]). These singularities are characterized by the 
fact that while the Hubble parameter and scale factor 
remain finite, all higher derivatives of the scale factor 
(a etc.) diverge as the cosmological singularity is ap- 
proached. These singularities may be obtained as the 
result of embedding of (3+l)-dimensional brane in the 
bulk and this is why these singularities will be hence- 
forth referred to as the "brane-like" singularities. We'll 
define the "brane-like" singularities in a following fash- 
ion: we'll say that singularity is of a "brane-like" type 
if at the instance of its occurrence both scale factor and 



density remain finite and nonzero, while all the higher 
order derivatives of scale factor ( starting with the second 
order) become altogether singular, i.e. a — > a s , p — > p s , 
< a s < oo, < p s < oo, d n a s /dt n — oo for n > 1. 

Evidently, the class of "brane-like" singularities in- 
cludes the singularities of Type II (with p s ^ 0) or 
SFS and BBtS. Moreover, BBS will also be of this type 
whenever we are talking about the models with the con- 
stant positive curvature, since at the singularity point 
p s = l/a 2 s . 

The physical nature of "brane-like" singularities emer- 
gence is quite clear: in the simple case with Z2 reflection 
symmetry and the identical cosmological constants on the 
two sides of the brane, the dynamical equation contains 
few additional terms. One of them is the square root of 
the sum of contributions of density (on the brane), ten- 
sion, cosmological constant and the "dark radiation" (the 
last one arises due to the projection of the bulk gravita- 
tional degrees of freedom onto the brane [1] ) . This sum is 
not positively defined and might become negative during 
the cosmological evolution. Thus, the solution of the cos- 
mological equations can't be continued beyond the point 
where this sum turns to zero and what we end up with 
at this point is nothing but a "brane-like" singularity. 
Since the existence of such singularities is natural in the 
brane physics, it won't be against the logic to assume 
that the appearance of " brane-like" singularities in usual 
Friedmann cosmology (SFS or BBtS) might be an evi- 
dence of validity of the brane hypothesis. Therefore it is 
interesting to consider "brane-like" singularities without 
a brane (i.e. in FLRW cosmology) to establish the par- 
ticular form of potential and the equation of state that 
will result in such singularities during cosmological dy- 
namics. Such potential and the equation of state may 
altogether be useful for answering the big cosmological 
question: Don't we really live on the brane? 

Furthermore, such singularities may actually result in 
very unusual models. In fact, let's consider the universe 
which contains a " brane-like" singularity. If the universe 
is filled with a scalar field while the Hubble parameter 
H(t s ) = H s and the scale factor a(t s ) = a s are finite at 
the singular point (H s < 00, a s < 00) then the value 
of the scalar field (f>(t s ) — cf> s might be finite as well. 
On the other hand, quantum corrections of higher order 
(in N-loops approximation) depend on the higher deriva- 
tives. If higher derivatives of scale factor diverge then 
this will also be the case for the scalar field. So one can 
expect that since all higher derivatives of scale factor and 
field alike diverge as the cosmological singularity is ap- 
proached, then the quantum effects will be dominating 
for t — > t s . This will be the case in spite of the fact that 
both density p s and scale factor a s will be finite and that 
p s might be small and a s - very large. 

It may seem that quantum corrections will be domi- 
nating because the pressure \p\ — > 00 as the cosmological 
singularity is approached. This is not the case for the sin- 
gularities of the IV type being "brane-like" by definition. 
Moreover, in Sec. VI we'll construct the singularities of 
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even mo re g eneral type that will violate the classifica- 
tions of [12| since p s and p s will be finite and nonzero 
and all higher derivatives will diverge. 

In this paper one constructs "brane-like" singulari- 
ties in Friedmann-Lemaitre-Robertson- Walker universe 
filled with the usual self-acting, minimally coupled scalar 
field or homogeneous tachyon field. In this cosmology 
we'll also construct the singularities (with finite scale fac- 
tor) where that Hubble variable vanishes and all higher 
derivatives of the scale factor diverge as the cosmologi- 
cal singularity is approached. That type of singularities is 
the generalization of the Big Break singularities and there 
will also be those of the "brane-like" type for the case of a 
constant positive curvature. We will calculate both self- 
acting potential V(<p) and tachyon potential V(T) that 
result in appearance of such singularities. Additionally, 
we'll present the equation of state for such models. More- 
over, the classification of singularities from the Ref. [l2[ 
will be complemented. 

This paper contains a discussion of a simple but use- 
ful method which allows one to construct exact solutions 
of the cosmological Friedmann equations filled with both 
self-acting, minimally coupled scalar field and a homoge- 
neous tachyon field. The method itself will be denoted 
as the method of linearization and it has been previously 
suggested in pj] (see also P3, 0, d). We'll give 
a brief discussion of this method in the next Section. 
In Sec. Ill we construct exact solutions with "brane- 
like" singularities. Then we consider field models (both 
tachyon and the usual minimally coupled scalar field) 
near the singularity (Sec. IV). Finally, in Sec. V we show 
that method of linearization allows one to obtain exact 
forms of potentials in tachyon models. In particular, we'll 
show that tachyon model that has been discussed in de- 
tail in Q is one of the simplest models in framework of 
the method of linearization. In Sec. VI we construct the 
new type of singularities which are some kind of gener- 
alization of the type IV singularities. The discussion is 
concluded in Sec. VII. 



II. THE METHOD OF LINEARIZATION 



Let us write the Friedmann equations as 



- =P-~, 2- = -(p + 3p). (2) 
a / er a 



The crucial point of this paper is the fact that a Fried- 
mann equations admits a linearizing substitution and can 
therefore be studied via the different powerful mathemat- 
ical methods which were specifically developed for the 
linear differential equations. This is the reason we call 
our approach the method of linearization [r| : 
Proposition. Let a = a(t) (with p = p(t), p = pit)) be 
a solution of ©. Then for the case k = the function 



ijj n = a n is the solution of the Schrodinger equation 

i — Un, 



(3) 



2 3n 
U n = n p - — {p + p) 



(4) 



with potential 
For example: 

Ui = - P+ ^ P , U 2 =p-3p, U 3 = ^(p-p) } 



TT 1 ( , 3 ^ TT 5/7 +3p 

and so on. 

Remark 3. If the universe is filled with scalar field 
whose Lagrangian is 



L= 2~ VW = K - y ' ( 5 ) 

then the expression (H|) will be 

U n = n{n - 3)K + n 2 V. 

In particular case n = 3 U3 — 9V(4>). This particular 
case has been extensively studied in 17], 18]. 
Remark 4- For small values of n 1 one gets U n ~ 
— 3n(p + p)/2; for example, if n = 0.01 then 



(0.0149p+0.015p) 0.015(p + p). 



Therefore one can use U n < to check whether the weak 
energy condition is violated plj . If, on the contrary, 
n 3> 1 then U n ~ n 2 p. 

Remark 5. If k = ±1 then the Proposition is valid only 
for the case n = 0, 1. 

As we have shown in jl9| , the representation of the 
Einstein-Friedmann equations as a second-order linear 
differential equation (J3|) allows for a usage of an arbi- 
trary (known) solution for construction of another, more 
general solution parameterized by a set of 3N constants, 
where N is an arbitrary natural number. The large num- 
ber of free parameters should prove itself useful for con- 
structing a theoretical model that agrees satisfactorily 
with the results of astronomical observations. In partic- 
ular, N = 3 solutions in the general case already exhibit 
inflationary regimes [l9j . Unlike the previously stud- 
ied two-parameter solutions (see [13], [3), t nese three- 
parameter solutions might describe an exit from inflation 
without any fine tuning of parameters as well as the sev- 
eral consecutive inflationary regimes. 

In the next Section we will show that the method of 
linearization is indeed an effective one for construction of 
a "brane-like" singularity. 



4 



III. "BRANE-LIKE" SINGULARITY 



one gets 



Assume 



U(t) = 



Hilt 



(ts - t) 



(6) 



with u 2 s = const >0, « = ±1, a > 0. For simplicity 
one has omitted the index n: U n — > U(t), ip n — > ip. For 
\t — t s \ >• 1 the potential U(t) — > so tp(t) ~ t and 
a(i) ~ t x / n . If n = 2 we have a universe filled with 
radiation (w — 1/3), and for n — 3/2 we have a dust 
universe with w = 0. 

Now let us consider the solution of the at t — > i s : 



fc=l 



where 



and 



= a y n = 



(2 -a) (a- 1) 



HUt 



(7) 



(8) 



(2-a)^k\Y[ k ~= [{k-m){2-a)-l]' 



with ci = 1. For a < 2 the series (JTJ) is convergent for 
any t (including t = t s ) since its radius is 

R = lim | — — |= +oo. 
If a > 2 then the first two terms of ([7]) are 

ijj = i/j s + 



(t s - ty 



so the function ip(t) will be singular at t — ¥ t s and we 
have either Big Rip (n > 0) or Big Crunch (n < 0) at 
t — t s . If a = 2 the general solution of the §Ji§ will be of 
the form 



ip = VtT^t (d(t s - t) L + c 2 (t s - ty L ) 



(10) 



with the arbitrary constants C\,i and L = v/l + 4/04^/2. 
For any C\ and C2 the solution (JTUJ) results in Big Rip 
or Big Crunch singularity as well. Finally, if a = 1 then 
the series will be 

</> = (!,- f) [1 + ^(*, - + f§(*. - <) 2 + 

so there is no cosmological singularity and a is finite. 

Thus, in framework of our investigation one must con- 
sider only < a < 1 and 1 < a < 2. Keeping only the 
two first terms in ([7]) and using 

tf = 4< '- = ^-(n-m\ (ID 

nip a nip 



1 (2 - a)(a - 1) 2 -a 



,2 



iJ = 



n(a — 1) 
0,2 



(12) 



a n(t s — i) c 



_ ul{n- 1) , _ ,2-a 

n{a — l) z 



Using (|12p one can show that the conditions ip s > and 
H > will hold if n > 0. So one has two cases: 

(i) If k = — 1 then < a < 1 and one gets BBS; 

(ii) if k = +1 then 1 < a < 2 and one gets BFS. 

To obtain SFS one should use another solution of (J3j) : 



ip = ip s - ntp s H s (t s — t) + 



(t s -tf- a + 



(2 -a)(l-a) 



C2*; (is - <) 



(2-a)(fe+l) 



fc=l 



+ C 2 fe+1 (t„ - t) 
k=0 



(2-a)(fc+l)+l 



(9) where if s = const > 0, tp s — al^ n — const > and 



(13) 



C2fe 



(k«2)* + V« 



(2-a)W(Hl)in^,K2-a)-l] 

(KU^) fc+1 n7J s ?/'s 
_ (2 - a)*+i (fc + l)inS N2 -") + !]' 

Using (fTTj) we get 

,2 



C2fc+1 



a \ ku s 
aj t ^ t n{t s -ty 



[n-l)Hi, 



so the solution (1131) contains the SFS at t — > t s for a < 1. 

At last, let's consider the ([7]). After differentiation one 
gets 



d m ip 
dt m 



00 m— 1 



= (-ir5> n [fc(2-«)-/](i s -<) fe 



k(2— a) — m 



k=l 1=0 

so we have a singularity in (|14[) for 



2-a 



(14) 
(15) 



Since < a < 1 then 



1 1 

2<2-^ <L 



Therefore for m > 1 the expression ([14]) diverge as the 
cosmological singularity is approached. Thus we have 
obtained the solution which contains a some kind of gen- 
eralization of the Big Break singularity (the scale factor 
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remains finite and the Hubble parameter vanishes as sin- 
gularity is approached) . In the case of positive curvature 
one has to choose n = 1 (cf. Remark 5) and we have a 
" brane-like" singularity (the density is finite and positive 
whereas all higher derivatives of the scale factor, starting 
out from the second one, diverge as the cosmological sin- 
gularity is approached). In the next Section we'll present 
a couple of models with the self-acting and minimally 
coupled scalar fields or with the homogeneous tachyon 
fields T — T(t) described by Sens or Born-Infeld type 
Lagrangians which result in such behavior. 

The similar investigation can be done for the ([13)1 . This 
singularity is characterized by the fact that Hubble pa- 
rameter remains finite instead of vanishing as the cosmo- 
logical singularity is approached. This solution describes 
the appearance of a "brane-like" singularity in the fiat 
universe. 



where w = p/p is defined by the fp~9|) . 

The expression for V holds iff w < 0. If w > 
one should introduce a new field theory based on a Born- 
Infcld type action with Lagrangian Q 



W(T)Vf 2 - 1, 



(23) 



so 



and 



W{T) 

Vf 2 - 1 ' 



p = L, 



ib 2 

T 2 = l + w, W=^—v'ir. 



(24) 



It is interesting that models ([20)) and (l23l) can be con- 
nected via the so called " transgression of the boundaries" 



IV. FIELD MODELS 

If the universe is filled with a self-acting and minimally 
coupled scalar field with Lagrangian ([5]) then the energy 
density and pressure are 



Now, using the potential ^ and the solution ([7} we 



get 



2n(a - If 
3nu 2 (t s -t) 2 



(25) 



therefore 



p = K + V, p = K-V, 



V=\{p-p), K=±(p+p). (16) 



Using ^ one can write 

tp 2 — ipxji (w + l)ip 
3nib 2 2n 2 tp 2 ' 



K = 



V 



mpib + (3 — n)ip (1 — w)ijj 2 



3n 2 i\) 2 



2n 2 %b 2 



where 



P i , 2n I i v>v> 

w = - = — 1 H 1 — 

p 3 1 i, 2 



(17) 



(18) 



(19) 



The second model is the universe filled with a homo- 
geneous tachyon field T = T{t) described by the Sens 
Lagrangian density [20]: 



L = -V(T)Jl-g 00 T 2 , 



(20) 



and in a flat Friedmann universe with metric ds 2 = dt 2 - 
a 2 (t)dr 2 we have density and pressure 



P = 



V{T) 



P 



(21) 



Therefore, for the case k — +1 we have w < (BFS) and 
w —} — oo as the cosmological singularity is approached. 
Using (fT7]) and (f22[) one concludes that this will be the 
case for the phantom (both usual (f> and tachyon T) fields 
only - the case which we leave out of consideration in this 
paper. 

For the case k = — 1 and < a < 1 one has a Big 
Break singularity with w — > +oo. For the model |[5} we 
get 



2u 2 



3n(t s - t) a 



>0, 



and 



2ui 



P 



3n(t s - t) a 



-> +oo, H ->• 0, 



as the cosmological singularity is approached. So all en- 
ergy conditions arc satisfied and H s — as it should be 
for a BBS. If k — +1 we have a "brane-like" singularity 
with 



Ps = 



(2-a) 2 (l-a) 2 



(26) 



At t — > t s the potential V = V(<j>) and field 4> = cf>(t) are 
given by expressions 



v{<j))= _^f3n(2-a) 



3n \ 8u s 



-a/(2-a) 



.)■ 



-2a/(2-a) 

(27) 



T 2 = 1 + w, 



V 




(22) 



2u s 
2 — a V 3n 



{ts-t) l - a/ \ 
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and V((f>) — »• — oo, (j> — > 4> s at t —> t s . 

For the solution (TTBI we have the same potential ([27)1 
and 



2nui 



instead of (j2"5")l . 

In the case of tachyon cosmology one should use the 
model ([23]). It easy to see that f 2 > and 



T(t) = T s T 2 



2n 1 



(ts - t) 



a/2 



Near the singularity, the correct equation of state for 
the solution ([7]) in case of a positive curvature has the 
form which looks similar to (|2T))) : 



p + 3p = 7 (p s - p) 



-|A| 



where 



2(1 



and < a < 1. 



(31) 



TACHYON POTENTIALS 



W{$) 



3n n(l — q 



with 



1/3 an, ,„ „ , 
^2V^l-t (r " Ts) - 

Since < a < 1 then 2/a — 3 > —1. For example, if 
a = 2/5 then 



IU(T) = 



54^5/2 



(T - T s f , 



and for the a — 2/7 



VF(T) 



21%/6u^ 
12500n 7 /2 



(T-T s 



Finally, let us consider the equation of state which re- 
sults in the potential (f27|) . It was shown in pj that the 
equation of state 



P + P = IP 



(29) 



results in dynamics which might be described by the self- 
acting potential is the form 



v(<f>) = 



= -2/(X-X) _ 7q-2A/(A-1) j 



(30) 



with 



Q 



a/y(A-i)(0-0.) 



where 7 > 0, A > 1. When 



we have 



Unfortunately, this expression is just formally equivalent 
to potential (|2"T|) . In fact, the second term in (f3"U|) (which 
is the dominant one at <fi — > <fi s ) is exactly (f2~T)) if a — 
2A/(2A - 1), therefore for A > 1 we get 1 < a < 2. It 
means that p s = 00 and we have a singularity of the 
III type which, in the case of general position, is not a 
"brane-like" one. 



The method of linearization proves to be extremely 
useful in finding the potentials of the exact solvable 
tachyon models. In particular, as we shall see, the 
tachyon model which was discussed in detail in Q is one 
of the simplest models in framework of the method of 
linearization. 

Let's consider Eq. ([3]) with potential U — 0. The solu- 
tion of this equation ift = Ct+C — > Ct by the translation 
t — > t — C'/C. In this case we get 



P 

w = — 



2n 

- 1 + T< 



v : 



2n- 3 
3nH 2 ' 



T 



w - 



1. 



so 



and 



±\/-3-(*-*.)+T„ 



V(T) 



2V9 - 6n 



n[2(T-T a )±y/Ent 8 ]' 



(32) 



If n = 3(1 + k)/2 with -1 < k < +1 (so < n < 3) 
and t s = then (1321) has exactly the form of one of the 
potentials from the first paper [7| (with To — > T s ). 

The case with U = p? = const > is a more inter- 
esting example . The solution of the (]3]) with Big Bang 
singularity at t = is ip = C sinh(/zi) . So 



w 



-1 



2n 



3 cosh p,t 
Using (|2"2"j) one gets 



P = 



p 2 (2n — 3 cosh 2 pt) 
2n 2 sinh 2 pt 



I 2n 

T(t) = ± J — arctan(sinh p,t) + T , 



and 



. . p 2 cosh at I n 

y W = 2 • u2 , V 9cosh ^"6n. 
3n z sinn pt 

Introducing A = p? /n 2 , k = 2n/3 — 1 we get 

V(T) = -4-v/l-(l + fc)cos 2 ^ 
sin 4 



(33) 
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where 



At f = f. 



£=-y/A(l + k)T. 

([55]) is the basic tachyon model of the paper Q . 

It is not difficult to construct many other integrable 
tachyon models using the simple, solvable potentials of 
d3j. Another fruitful way of doing it lies in a use of 
the Darboux transformation to those initial potentials 
(U = 0, U = p 2 ). This, however, is out of scope of the 
paper. 



VI. GENERALIZED TYPE IV SINGULARITIES 

The singularities of type IV (according to the classi- 
fication of Ref. [T3|) have the following behavior: for 
t —> t s , a —} a s , p — > 0, \p\ —> and the higher derivatives 
of H diverge (0 < a s < oo). In this section we present a 
new type of singularities: 

Generalized type IV: For t — > t s , a — > a s , p — > p Sl 
p — > p s and higher derivatives (starting out from the 
third one) of H diverge and 0<a s <oo,0</9s<oo, 
< \p s \ < oo. 

Let's put n = 3, ip = a 3 and (in parametric form) 

ip = A + ^j- (4 cos 7] — 2 cos 2 Ti — cos 4 r\ — 4 log(l + cos 77)) , 



t = t. 



- (log 



1] 

tan- 



cost? 



(34) 

where A, B, k are constants, 0<T7<7r;?7 = corre- 
sponds to t = — oo, 77 = 7r/2 to t = t s (singularity) and 
rj = 7r to t = +00. After the differentiation we get (a 
dot denotes the derivative with respect to cosmic time t 
rather than to parameter 77) 

%Ij = k 2 B (l - cos 3 77) , 



■ip = 3k 3 B sin 2 77, 



6Bk sin 77 
cos 77 



tp = A, ijj = k 2 B 7 ip = 3k 3 B 

and, starting out from the third one, all higher derivatives 
diverge. At t = +00 

ip = sign(K-B) x 00, ip — 2k 2 B, tp = 0. 

Therefore at t = t s 



K 4 B 2 
9A 2 



Ps 



Ps 



Bk 3 (kB - 6 A) 
9A 2 : 



1 



6A 
~kB' 



Thus we have a generalization of a type IV singularity at 
t s , where the density and pressure are finite and nonzero 
whereas all higher derivatives of H diverge. 

It is convenient to introduce a new parameter s: 

QA 



so 



4k 2 

Ps = —j 



Ps 



4k 2 (1- s) 



1 - s. 



If 4/3 < s < 2 then -1 < w s < -1/3; if s = 2 then 
w s = — 1; if s > 2 then w g < — 1. To obtain the initial 
Big Bang singularity at t — U, —00 < ti < t s one should 
put 



s < Si = 6 log 2 - - - 2.659, 

or 



w s > Wi = —1.659. 

In the initial Big Bang singularity, the barotropic index 
77; = +00, on the other hand, w{j]) is the monotonously 
decreasing function. These properties result in a follow- 
ing conclusion: if 

4 

g < s < Si, 

then after Big Bang the model will go through the 
usual expansion with damping, but starting out from 
some moment it will experience an accelerated expansion 
up to a future generalized type IV singularity. 



6k 5 B sin 2 77(cos 2 7/ + 1) 

V = 4 , 

cos 4 77 

and so on. 

For kB > the function (pH)) is the monotonously 
increasing one for < 77 < tt and ip{^) — +00 (i.e. at 
t = +00). Thus at t = —00 

ip = A + Bk Q -log2J , t/> = 0, tA = 0. 



VII. CONCLUSION 

In this paper we have discussed a simple method of 
construction of exact solutions of the Friedmann equa- 
tions with finite scale factor singularities. Despite sim- 
plicity, the method allows for acquirement of solutions 
characterized by the extremely interesting properties. 

The main results of this work are: 

(i) we have obtained a new type of finite-time, future 
singularities which seem to be most similar to the type IV 



8 



of [12[ but are different nevertheless as they have nonzero 
pressure and density at the singular point; 

(ii) we have obtained a new type of finite-time, fu- 
ture quasi-singularities being rather similar to the w- 
singularities (which are quasi-singularities too) but hav- 
ing nonzero pressure at a singular point; 

(iii) we have shown that "brane-like" singularities can 
occur in a common Friedmann cosmology with potential 
(l27|) and equation of state (|3Tj) (near of singularity) as 
well as in a tachyon cosmology with potential (|28|) ; 

(iv) we have obtained the generalized Big Break singu- 
larities not only for the universe filled with tachyons but 
also a usual minimally coupled scalar field; 

(v) we have shown that basic tachyon model which was 
discussed in detail in Q is one of the simplest models in 
framework of the method of linearization. 
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